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Modern-day 'testing' of (perturbative) QCD is as much about pushing 
the boundaries of its applicability as about the verification that QCD is the 
correct theory of hadronic physics. This talk gives a brief discussion of a 
small selection of topics: factorisation and jets in diffraction, power correc- 
tions and event shapes, the apparent excess of 6-production in a variety of 
experiments, and the matching of event generators and NLO calculations. 

PACS numbers: 12.38.-t, 12.38.Aw 
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1. Introduction 

The testing of QCD is a subject that 
many would consider to be well into matu- 
rity. The simplest test is perhaps that a s 
values measured in different processes and 
at different scales should all be consistent. 
It suffices to take a look at compilations 
by the PDG [2] or Bethke [3] to see that 
this condition is satisfied for a range of ob- 
servables, to within the current theoretical 
and experimental precision, namely a few 
percent. There exist many other poten- 
tially more discriminatory tests, examples 
explicit measurements of the QCD colour 
factors [4] or the running of the 6-quark 
mass [5] — and there too one finds a sys- 
tematic and excellent agreement with the 
QCD predictions. A significant amount of 
the data comes from HERA experiments, 
and to illustrate this, figure 1 shows a com- 
pilation of a subset of the results on a s , as compiled by ZEUS [1]. 
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H — Experimental errors 
— Theoreticol error 

ZEUS (prel. ) 96-97 Incl. jets 

ZEUS NLO-QCD fit (prel.) 2001 

ZEUS (prel.) 96-97 Subjets 

ZEUS (prel.) 95-97 Jet shope 

ZEUS dijets 

(Phys. Lett. B507 (2001) 70) 
HI Incl. jets 

(Eur. Phys. J.C19 (2001) 289) 

H1 NLO-QCD fit 
(hep-ex/0012053) 

PDG 2000 

(Eur. Phys. J. C15 (2000) 1) 



(J. Phys. G26 (2000) R27) 
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Fig. 1. A compilation of HERA 
a s measurements, taken from [ll. 
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In the space available however, it would be impossible to give a critical 
and detailed discussion of the range of different observables that are used 
to verify that QCD is 'correct'. Rather let us start from the premise that, 
in light of the large body of data supporting it, QCD is the right theory of 
hadronic physics, and consider what then is meant by 'testing QCD'. 

One large body of activity is centred around constraining QCD. This 
includes such diverse activities as measuring fundamental (for the time be- 
ing) unknowns such as the strong coupling and the quark masses; measuring 
quantities such as structure functions and fragmentation functions, which 
though formally predictable by the theory are beyond the scope of the tools 
currently at our disposal (perturbation theory, lattice methods); and the 
understanding, improvement and verification of the accuracy of QCD pre- 
dictions, through NNLO calculations, resummations and projects such as 
the matching of fixed-order calculations with event-generators. One of the 
major purposes of such work is to provide a reliable 'reference' for the inputs 
and backgrounds in searches for new physics. 

A complementary approach to testing QCD is more about exploring the 
less well understood aspects of the theory, for example trying to develop an 
understanding of non-perturbative phenomena such as hadronisation and 
diffraction, or the separation of perturbative and non-perturbative aspects 
of problems such as heavy-quark decays; pushing the theory to new limits 
as is done at small- x and in studies of saturation; or even the search for and 
study of qualitatively new phenomena and phases of QCD, be they within 
immediate reach of experiments (the quark-gluon plasma, instantons) or not 
(colour superconductors)! 

Of course these two branches of activity are far from being completely 
separated: it would in many cases be impossible to study the less well 
understood aspects of QCD without the solid knowledge that we have of 
its more 'traditional' aspects — and it is the exploration of novel aspects of 
QCD that will provide the 'references' of the future. 

The scope of this talk is restricted to tests involving final states. Final 
states tend to be highly discriminatory as well as complementary to more 
inclusive measurements. We shall consider two examples where our under- 
standing of QCD has seen vast progress over the past years, taking us from 
a purely 'exploratory' stage almost to the 'reference' stage: the question of 
jets and factorisation in diffraction (section 2); and that of hadronisation 
corrections in event shapes (section 3) . We will then consider two questions 
that are more directly related to the 'reference' stage: the topical issue of 
the excess of 6-quark production seen in a range of experiments (section 4); 
and then the problem of providing Monte Carlo event generators that are 
correct to NLO accuracy, which while currently only in its infancy is a 
subject whose practical importance warrants an awareness of progress and 
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pitfalls. 

For reasons of lack of space, many active and interesting areas will not be 
covered in this talk, among them small- x physics, progress in next-to-next- 
to-leading order calculations, questions related to prompt photons, the topic 
of generalised parton distributions and deeply-virtual Compton scattering, 
hints (or not) of instantons, a range of measurements involving polarisation 
and so on. Many of these subjects are widely discussed in other contribu- 
tions to both the plenary and parallel sessions of this conference, to which 
the reader is referred for more details. 

2. Jets in diffraction and factorisation 

Factorisation, for problems explicitly involving initial or final state hadrons, 
is the statement that to leading twist, predictions for observables can be 
written as a convolution of one or more non-perturbative but universal 
functions (typically structure or fragmentation functions) with some per- 
turbatively calculable coefficient function. 

While factorisation has long been es- 
tablished in inclusive processes [7] it has 
been realised in the past few years [8] 
that it should also hold in more exclusive 
cases — in particular for diffraction, in 
terms of diffractive parton distributions 
fa)p ( x ' x iP' A* 2 ' which can be interpreted 
loosely as being related to the probabil- 
ity of finding a parton a at scale /i 2 with 

longitudinal momentum fraction x, inside Fig. 2. Illustration of diffractive 
a diffractively scattered proton p, which kinematics. Figure taken from [6]. 
in the scattering exchanges a squared mo- 
mentum t and loses a longitudinal momentum fraction . These kinematic 
variables are illustrated in fig. 2. 

The dependence of the diffractive parton distributions on so many vari- 
ables means that without a large kinematical range (separately in x, x-p and 
Q 2 , while perhaps integrating over t) it is a priori difficult to thoroughly 
test diffractive factorisation. An interesting simplifying assumption is that 
of Regge factorisation, where one writes [9] 

/ a d ;J(x,x P ,/i 2 ,t) = \f5 p {t)\ 2 x^ a{t) f a/w {x/x w ^ 2 ,t) (1) 

the interpretation of diffraction being due to (uncut) pomeron exchange 
(first two factors), with the virtual photon probing the parton distribution 
of the pomeron (last factor). 

As yet no formal justification exists for this extra Regge factorisation. 
Furthermore given that diffraction is arguably related to saturation and 
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high parton densities (assuming the AGK cutting rules [10]) one could even 
question the validity of arguments for general diffractive factorisation, which 
rely on parton densities being low (as does normal inclusive factorisation). 

The experimental study of factorisation in diffraction relied until re- 
cently exclusively on inclusive Frf measurements. This was somewhat un- 
satisfactory because of the wide range of alternative models able to repro- 
duce the data and even the existence of significantly different forms for the 
ja/p (x I xp , /i 2 , t) which gave a satisfactory description of the data within 
the Regge factorisation picture. 



H1 Diffractive Dijets 
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Fig. 3. Comparisons of HI diffractive dijet cross sections with predictions obtained 
using the assumption of Regge factorisation [11]. 



However diffractive factorisation allows one to predict not only inclusive 
cross sections but also jet cross sections. Results in the Regge factorisation 
framework are compared to data in figure 3 (taken from [11]), showing 
remarkable agreement between the data and the predictions (based on one 
of the pomeron PDF fits obtained from F^). On the other hand, when 
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one considers certain other models that work well for the disagreement 
is dramatic, as for example is shown with the soft colour neutralisation 
models [12, 13] in figure 4. 

Despite this apparently strong confir- 
mation of diffractive factorisation, a word 
of warning is perhaps needed. Firstly there 
exist other models which have not been 
ruled out (for example the dipole model 
[14]). In these cases it would be of in- 
terest to establish whether these models 
can be expressed in a way which satisfies 
some effective kind of factorisation. 

Other important provisos are that a 
diffractive PDF fit based on more recent z (i<= ls ) 

F$ data has a lower gluon distribution lp 
and so leads to diffractive dijet predic- 
tions which are a bit lower than the data, 
though still compatible to within experi- 
mental and theoretical uncertainties [15]. 
And secondly that the predictions them- 
selves are based on the Rapgap event gen- 
erator [16] which incorporates only leading order dijet production. It would 
be of interest (and assuming that the results depend little on the treat- 
ment of the 'pomeron remnant,' technically not at all difficult) to calculate 
diffractive dijet production to NLO with programs such as Disent [17] or 
Disaster++ [18], using event generators only for the modelling of hadroni- 
sation correction, as is done in inclusive jet studies. 



Fig. 4. A comparison between 
diffractive dijet data and results 
from the soft-colour interaction 
(SCI) [12] and semiclassical [13] 
models. Figure adapted from [11]. 



3. Hadronisation 

Another subject that has seen considerable experimental and theoretical 
progress recent years is that of hadronisation. Even at the relatively high 
scattering energies involved at LEP and the Tevatron, for many final state 
observables non-perturbative contributions associated with hadronisation 
are of the same order of magnitude as next-to-leading order perturbativc 
contributions and cannot be neglected. With the advent of NNLO cal- 
culations in the foreseeable future the need for a good understanding of 
hadronisation becomes ever more important. 

Until a few years ago, the only way of estimating hadronisation correc- 
tions in final-state measurements was by comparing the parton and hadron 
levels of Monte Carlo event generators. Such a procedure suffers from a 
number of drawbacks. In particular the separation between perturbative 
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and non-perturbative contributions is ill-defined: for example event genera- 
tors adopt a prescription for the parton level based on a cutoff; on the other 
hand, in fixed-order perturbative calculations no cutoff is present, and the 
perturbative integrals are naively extended into the non-perturbative region 
— furthermore the 'illegally-perturbative' contribution associated with this 
region differs order by order (and depends also on the renormalisation scale) . 

Additionally, hadronisation corrections obtained from event generators 
suffer from a lack of transparency: the hadronisation models are generally 
quite sophisticated, involving many parameters, and the relation between 
these parameters and the hadronisation corrections is rarely straightforward. 

In the mid 1990's a number of groups started examining approaches for 
estimating hadronisation corrections based on the perturbative estimates of 
observables' sensitivity to the infrared. This leads to predictions of non- 
perturbative corrections which are suppressed by powers of 1/Q relative 
to the perturbative contribution (for a review see [19]). One of the most 
successful applications of these ideas has been to event shapes, for which 
(in the formalism of Dokshitzer and Webber [20] ) 

2Cf in 



(V N p) = (Vpt) + c v V . 



V = 



TT 
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where cy is a perturbatively calculable 0.7 
observable-dependent coefficient and V gov- 
erns the size of the power correction. The 
quantity ao(fj,i), which can be interpreted 
as the mean value of an infrared finite ef- 
fective coupling in the infrared (up to an : 
infrared matching scale convention- 
ally chosen to be 2 GeV), is hypothesised 
to be universal. The terms in powers of 
a s are subtractions of pieces already in- 
cluded in the perturbative prediction for 
the observable. 

It is interesting to see the progress 
that has been made in our understand- 
ing of these effects. The first predictions 
for the cv coefficients were based on cal- 
culations involving the Born configuration plus a single 'massive' (virtual) 
gluon. Fitting ao and a s to data for mean values of e + e~ event-shapes, using 
the original predictions for the cy, leads to the results shown in figure 5. 

At the time of the original predictions, however, much of the data used 
to generate fig. 5 was not yet in existence (which is perhaps fortunate - 
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Fig. 5. Results of fits to e + e~ mean 
event shapes using original, 'naive' 
calculations for Cy. 
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had fig. 5 been around in 1995, the field of 1/Q hadronisation corrections 
might not have made it past early childhood). Rather, various theoretical 
objections (e.g. [21]) and the gradual appearance of new data, especially for 
the broadenings, forced people to refine their ideas. 

Among the developments was the re- 
alisation that to control the normalisation 
of the cy it is necessary to take into ac- 
count the decay of the massive, virtual, 
gluon (the reason for the two thrust re- 
sults in fig. 5 was the existence of two 
different conventions for dealing with the 
undecayed massive gluon) [22]. It was 
also realised that it is insufficient to con- 
sider a lone 'non-perturbative' gluon, but 
rather that such a gluon must be taken 
in the context of the full structure of soft 
and collinear perturbative gluon radiation 
[23] . Another discovery was that hadron- 
masses can be associated with universal- 
ity breaking 1/Q power corrections in cer- 
tain definitions of observables [24] and when 

testing the universality picture all observables should be measured in an ap- 
propriate common 'hadron-mass' scheme. 

Results incorporating these theoretical developments are shown in fig- 
ure 6. As well as e + e~ mean event shapes we also include recent results 
using resummed DIS event shapes [25], fitted to HI distributions [26]. The 
agreement between observables, even in different processes, is remarkable, 
especially compared to fig. 5, and a strong confirmation of the universality 
hypothesis. 1 

This is not to say that the field has reached maturity. In the above fits 
the approximation has been made that non-perturbative corrections just 
shift the perturbative distribution [29], however there exists a considerable 
amount of recent work which examines the problem with the more sophis- 
ticated 'shape-functions' approach [30] in particular in the context of the 
Dressed Gluon Exponentiation approximation [31]. An important point also 
is that all the detailed experimental tests so far are for 2-jet event shapes, 
where there exists a solid theoretical justification based on the Feynman 
tube model [32], i.e. longitudinal boost invariance. It will be of interest to 
see what happens in multi-jet tests of 1/Q hadronisation corrections where 

1 It should be noted that results for certain e + e~ distributions [27] and DIS means 
[26, 28] are not quite as consistent. Though this remains to be understood, it may in 
part be associated with the particular fit ranges that are used. 
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one introduces both non-trivial geometry and the presence of gluons in the 
Born configuration [33] . Finally we note the provocative analysis by the Del- 
phi collaboration [34] where they show that a renormalisation-group based 
fit prefers an absence of hadronisation corrections, at least for mean values 
of event shapes, as well as leading to highly consistent values for a s across 
a range of event-shapes. 



4. Heavy quark (b) production 

For light quarks (and gluon) it is impossible to make purely perturbative 
predictions of their multiplicity or of their fragmentation functions because 
of soft and collinear divergences. For heavy quarks however, these diver- 
gences are cut off by the quark mass itself, opening the way to a range of 
perturbative predictions and corresponding tests of QCD. 
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Fig. 7. Left: &-quark p t distribution at the Tevatron [35]; upper right: summary 
of open b cross sections in -fp, DIS and 77 collisions, normalised to theoretical 
expectations (figure taken from [36]); lower right: ratio of experiment to theory for 
the charm p t distribution at HERA (taken from [36]). 



It is therefore particularly embarrassing that there should be a signifi- 
cant discrepancy in most experiments (but not all, e.g. [37]) where the QCD 
bottom production cross section has been measured. The situation is shown 
in figure 7 for Tevatron, HERA and LEP results, illustrating the systematic 
excess of a factor of three between measurements and NLO calculations. 
To add to the puzzle, the agreement for charm production (which if any- 
thing should be worse described because of the smaller mass) is considerably 
better across a range of experiments (see e.g. the lower-right plot of fig. 7). 
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Aside from the intrinsic interest of having a good understanding of b- 
production in QCD, one should keep in mind that 6-quarks are widely relied 
upon as signals of Higgs production and in searches for physics beyond the 
standard model, so one needs to have confidence in predictions of the QCD 
background. 

We shall discuss a couple of explanations that have been proposed for the 
excess at the Tevatron (the excesses in other experiments are more recent 
and have yet to be addressed in the same detail). Indeed, one hypothesis 
is precisely that we are seeing a signal of light(ish) gluino production. An- 
other is that bottom fragmentation effects have been incorrectly accounted 
for. A third explanation, discussed in detail in another of the opening ple- 
nary talks [38] is associated with unintegrated kt distributions and small- x 
resummations. 



f.l. The SUSY hypothesis 




In [39] it has been argued that a pos- 
sible explanation of the Tevatron 6-quark 
excess is the production of a pair of light 
gluinos with a mass of order 14 GeV which 
then decay to sbottoms (~ 3.5 GeV) and 
bottoms, as in fig. 8. The mixing angles 
are chosen such that the sbottom decou- 
ples from the Z at LEP, accounting for 
its non-observation there. 

At moderate and larger pt , the contri- 
bution from this process is about as large 
as that from NLO QCD and so it brings 
the overall production rate into agreement 
with the data. 

There are a number of other conse- 
quences of such a scenario: one is the 
production of like-sign b quarks (as in the 
Feynman graph of fig. 8), which could in 
principle be observed at the Tevatron, al- 
though it would need to be disentangled 
from Bq-Bq mixing. Another is that the 
running of a s would be modified signifi- 
cantly above the gluino mass, leading to 
an increase of about 0.007 in the running 

to Mz of low Q measurements of a s . This seems to be neither favoured nor 
totally excluded by current a s measurements. 
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Though they have not provided a detailed analysis, the authors of 
also consider the implications for HERA. There it seems that the enhance- 
ment of the 6-production rates is too small to explain the data (because of 
the suppression due to the gluino mass). 

4-2. The fragmentation explanation 

In any situation where one sees a significant discrepancy from QCD 
expectations it is worth reexamining the elements that have gone into the 
theoretical calculation. Various groups have considered issues related to 6 
fragmentation and found significant effects, which could be of relevance to 
the Tevatron results (see for example [40]). However a recent article by 
Cacciari and Nason [41] is particularly interesting in that it makes use of 
the full range of available theoretical tools to carry out a unified analysis all 
the way from the e + e~ data, used to constrain the 6-quark fragmentation 
function, through to expectations for the Tevatron. It raises a number of 
important points along the way. 2 

To be able to follow their analysis it is worth recalling how one calculates 
expectations for processes involving heavy quarks. The cross section for 
producing a 6-quark with a given p t (or even integrated over all pt) is finite, 
unlike that for a light quark. This is because the quark mass regulates (cuts- 
off) the infrared collinear and soft divergences which lead to infinities for 
massless quark production. But infrared finiteness does not mean infrared 
insensitivity and to obtain a l?-meson pt distribution from a 6-quark pt 
distribution, one needs to convolute with a fragmentation function, 

^- = f dp t dz ^- D(z) 5{p t - zp t ) . (3) 

apt J apt 

^^s-' ^-v^" fragmentation 

measured, e.g. _Bo PT QCD, b quark 

The details of the infrared finiteness of the 6-quark production are such that 
(zD(z)) is 1 — O (A/ nib), where the origin of the A/m^ piece is closely related 
to that of the A/Q power corrections discussed in the previous section [42]. 

There are various well-known points to bear in mind about fragmen- 
tation functions. Firstly, in close analogy to the hadronisation corrections 
discussed earlier (and of course structure functions), the exact form for the 
fragmentation function will depend on the perturbative order at which we 
define eq. (3). Secondly, while for pt ~ m& we are free to use fixed order 
(FO) perturbative predictions, for pt 3> there are large logarithmically 
enhanced terms, which need to be resummed. The technology for doing this 
currently exists to next-to-leading logarithmic (NLL) order. In the interme- 



2 The reader is referred to their article for full references to the 'ingredients' used at 
different stages of the analysis. 
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diate region pt > nib the two approaches can be combined to give FONLL 
predictions [43, 44] (strictly this can be used even for p t 3> nib). 

Having established these points we can consider what has been done by 
Cacciari and Nason [41]. Firstly they discuss moments of the fragmentation 
function (z N ~ l D{z)). This is because for a steeply falling perturbative pt 
distribution in eq. (3), ~ ^/Pt \ after integrating out the 5-function to 
give pt = pt/z, one obtains the result 



da 
dpt 



D(z)) 



B— meson 



da 
dpt 



(4) 



b— quark 




solid: N=Z fit 
dashed: e = 0.006 
dotted: e = 0.002 
dot-dashed: Pure NLL pQCD 
, , , , I , L 



where for the Tevatron N ~ 5. 

The cleanest place to constrain b i 
fragmentation is in e + e~ collisions. Fig- 
ure 9 shows moments of the momen- 
tum fraction (with respect to Q/2) 
carried by i?-mesons as measured by *'g 
Aleph [45] . The (magenta) dot-dashed 
curve shows the purely perturbative 
NLL prediction, which is clearly above t 
the data. The dashed curve shows 
what happens when one includes the 1 5 10 15 

convolution with an e = 0.006 Peter- 
son fragmentation function [46] . Why 

this particular function? Simply be- Fi S- 9 - Moments of the momentum 
cause it is the one included in certain fraction carried by B-mesons in e+e , 
Monte Carlo event generators and used com Pared to NLL predictions with and 
widely by experimental collaborations without fragmentation functions [41] . 
that have compared measured and theoretical pt distributions. The data 
point for the N = 5 moment is 50% higher than the theoretical expectation 
with this fragmentation function. 

Of course we don't expect agreement: the e = 0.006 Peterson is widely 
used in Monte Carlos where one has only leading-logs. But we are interested 
in NLL calculations and the fragmentation function needs to be refitted. 
The authors of [41] take the functional form of [47], fitted to the N = 2 
moment, to give the solid curve. 

The next step in the Cacciari and Nason analysis should simply have 
been to take the FONLL calculation of bottom production at the Tevatron 
[43], convolute with their new fragmentation function and then compare 
to data. This however turns out to be impossible for most of the data, 
because it has already been deconvoluted to 'parton-levek (in some cases 
with the e = 0.006 Peterson fragmentation function). So they are only able 
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Fig. 10. Left: result for B-meson production [48] compared to the FONLL pre- 
diction with the l N — 2' fragmentation function [41]. Right: results for 6-jet 
production [50] compared to the NLO predictions [51]. 



to compare with the recent CDF data [48] for l?-mesons, shown in the left- 
hand plot of figure 10. The dashed curve is the central result, while the solid 
ones are those obtained when varying the factorisation and renormalisation 
scales by a factor of two. 3 The dotted curve shows the results that would 
have been obtained with the Peterson fragmentation function. Predictions 
with FO (generally used in previous comparisons) rather than FONLL would 
have have been 20% lower still. 

Another interesting approach to the problem is to eliminate the frag- 
mentation aspects altogether, which can be achieved by looking at the Et 
distribution of 6-jets, without specifically looking at the b momentum [51]. 
This has been examined by the DO collaboration [50] and the comparison to 
NLO predictions is shown in the right-hand plot of figure 10. Though in a 
slightly different Et range, the relation between theory and data is similar to 
that in the Cacciari-Nason approach for l?-mesons: there is a slight excess 
in the data but not significant compared to the uncertainties. A minor point 
to note in the study of 6-jets is that there are contributions a" ln 2n_1 Et/rrib 
from soft and collinear logs in the multiplicity of gluons which can then 
branch collinearly to bb pairs [52]. At very large Et these terms would need 
to be resummed. 

So overall, once one has a proper theoretical treatment, including both 
an appropriate fragmentation function and, where relevant, an FONLL 
perturbative calculation, it is probably fair to say that the excess of b- 
production at the Tevatron is not sufficiently significant to be worrisome 

3 A point worth keeping in mind [49] is that the central scale choice /i = \Jp\ + m% 
is not universally accepted as being optimal — indeed for p t > mi,, a scale choice of 
fi — pt is equally justifiable, and would have a non-negligible effect on the predictions. 
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(or evidence for supersymmetry). 

At some of the other experiments where an excess of 6-production is 
observed a number of the same issues arise, in particular relative to the use 
of the e = 0.006 Peterson fragmentation function and the presentation of 
results at parton level rather than hadron level. However fragmentation is 
less likely to be able to explain the discrepancies, because of the lower p t 
range. 



5. Event generators at NLO 

The problem of matching event generators with fixed order calculations 
is one of the most theoretically active areas of QCD currently, and consid- 
erable progress has been made in the past couple of years. This class of 
problems is both of intrinsic theoretical interest in that it requires a deep 
understanding of the structure of divergences in QCD and of phenomenolog- 
ical importance because of the need for accurate and reliable Monte Carlo 
predictions at current and future colliders. 

Two main directions are being followed: one is the matching of event- 
generators with leading-order calculations of ra-jet production (where n may 
be relatively high), which is of particular importance for correctly estimating 
backgrounds for new-particle searches involving cascades of decays with 
many resulting jets. For a discussion of this subject we refer the reader 
to the contributions to the parallel sessions [53]. 

The second direction, still in its infancy, is the matching of event gen- 
erators with next-to-leading order calculations (currently restricted to low 
numbers of jets), which is necessary for a variety of purposes, among them 
the inclusion of correct rate estimates together with consistent final states, 
for processes with large NLO corrections to the Born cross sections (e.g. K 
factors in pp and jp collisions, boson-gluon fusion at small- x in DIS). 

While there have been a number of proposals concerning NLO matching, 
many of them remain at a somewhat abstract level. We shall here concen- 
trate on two approaches that have reached the implementational stage. As 
a first step, it is useful to recall why it is non-trivial to implement NLO 
corrections in an event generator. Let us use the toy model introduced by 
Frixione and Webber [54], involving the emission only of 'photons' (simpli- 
fied, whose only degree of freedom will be their energy) from (say) a quark 
whose initial energy is taken to be 1. For a system which has radiated n 
photons we write a given observable as 0(E q , E yi , . . . , E Jn ). So for exam- 
ple at the Born level, the observable has value 0(1). At NLO we have to 
integrate over the momentum of an emitted photon, giving the following 
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contribution to the mean value of the observable: 



f dx 

a / — R(x) 0(1 — x,x) , 



(5) 



where R(x) is a function associated with the real matrix element for one- 
photon emission. There will also be NLO virtual corrections and their 
contribution will be 



where V{x) is related to the matrix element for virtual corrections. 

The structure of dx/x divergences is typical of field theory. Finite- 
ness of the overall cross section implies that for x — > 0, R{x) = V{x). 
This means that for an infrared safe observable (i.e. one that satisfies 
lim x _»oO(l — x,x) = 0(1)), the O (a) contribution to the mean value of 
the observable is also finite. However any straightforward attempt to imple- 
ment eqs. (5) and (6) directly into an event generator will lead to problems 
because of the poor convergence properties of the cancellation between di- 
vergent positively and negatively weighted events corresponding to the real 
and virtual pieces respectively. So a significant part of the literature on 
matching NLO calculations with event generators has addressed question 
of how to recast these divergent integrals in a form which is practical for 
use in an event generator (which must have good convergence properties, 
especially if each event is subsequently going to be run through a detector 
simulation). The second part of the problem is to ensure that the normal 
Monte Carlo event generation (parton showering, hadronisation, etc.) can 
be interfaced with the NLO event generation in a consistent manner. 

One approach that has reached the implementational stage could be 
called a 'patching together' of NLO and MC. It was originally proposed in 
[55] and recently further developed in [56] and extended in [57]. There one 
chooses a cutoff x zcro on the virtual corrections such that the sum of Born 
and virtual corrections gives zero: 



It is legitimate to sum these two contribution because they have the same 
(Born) final state. Then for each event, a real emission of energy x is 
generated with the distribution dx/xR(x) and with the same cutoff as on 
the virtuals. The NLO total cross section is guaranteed to be correct by 
construction: 




(6) 
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The next step in the event generation is to take an arbitrary separation 
parameter x sep , satisfying x zcro < x scp < 1. For x > x sep the NLO emission 
is considered hard and kept (with ideally the generation of normal Monte 
Carlo showering below scale x, as in the implementation of [57]). For x < 
x sep the NLO emission is thrown away and normal parton showing is allowed 
below scale x. 4 

Among the advantages are that the events all have positive and uniform 
weights. And while the computation of x ZCTO is non-trivial, the method 
requires relatively little understanding of the internals of the event generator 
(which are often poorly documented and rather complicated). However the 
presence of the separation parameter x sep is in principle problematic: there 
can be discontinuities in distributions at x sep , certain quantities (for example 
the probability for a quark to have radiated an amount of energy less than 
some x r which is below x sep ) will not quite be correct to NLO and above 
^sc P potentially large logarithms of x sep are being neglected. These last two 
points mean that for each new observable that one studies with the Monte 
Carlo program, one should carry out an analysis of the x se p dependence 
(varying it over a considerable range, not just a factor of two as is sometimes 
currently done). 

A rather different approach (which we refer to as 'merging') has been 
developed by Frixione and Webber in [54]. 5 They specify a number of 
conditions that must be satisfied by a Monte Carlo at NLO (MC@NLO): 
i) all observables should be correct to NLO; ii) soft emissions should be 
treated as in a normal event generator and hard emissions as in an NLO 
calculation; iii) the matching between the hard and soft regions should 
be smooth. Their approach exploits the fact that Monte Carlo programs 
already contain effective real and virtual NLO corrections, 

±a^M(x) for ^ ual . (9) 

Because Monte Carlo programs are designed to correctly reproduce the 
structure of soft and collinear divergences, M(x) has 6 the property that 
for x — ► 0, M(x) = R(x) = V(x), i.e. the divergent part of the NLO cor- 
rections is already included in the event generator. This can be exploited 
when adjusting the Monte Carlo to be correct to NLO, because the regions 

4 For simplicity, many important but sometimes tricky technical details have been left 
out. This will also be the case for the merging procedure discussed lower down. 

5 A number of aspects of the work of Collins and collaborations [58] may actually 
be equivalent, though presented in a rather different framework. Related issues are 
discussed also in [59]. 

6 Or rather, 'should have.' In practice the divergence structure of large-angle soft-gluon 
emission is not always properly treated in event generators, which leads to some extra 
complications in the MC@NLO approach. 
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that need adjusting are the hard regions, but not the (soft) divergent re- 
gions. Specifically the method introduced in [54] can be summarised by the 
formula 

/dx 
— (V(x)-M(x)) 

/dx 
(10) 

-^MC,Born is to be read 'interface to Monte Carlo.' It means that one should 
generate a Monte Carlo event starting from the Born configuration (or from 
the Born configuration plus a photon in the case of iMC,Bom+x)- Since at 
the Born level, lMC,Bom already contains effective real and virtual correc- 
tions which go as ±aM(x)/x, when evaluating the NLO corrections to the 
MC, these pieces should be subtracted from the full NLO matrix elements. 
Because M(x) and R{x) (or V(x)) have the same x — ► limit, the real and 
virtual integrals are now individually finite and well-behaved, which means 
that the Monte Carlo only needs only a small, O (a), correction in order for 
it to be correct to NLO. 

Illustrative results from this ap- 
proach are shown in figure 11 for the 
transverse momentum distribution of 
a W + W~ pair in hadron-hadron col- 
lisions. In the low transverse momen- 
tum region (which requires resumma- 
tion — the pure NLO calculation breaks 
down) MC@NLO clearly coincides with 
the Herwig results, while at high trans- 
verse momentum it agrees perfectly 
with the NLO calculation (default Her- 
wig is far too low). 

So this procedure has several ad- Fig. 11. Transverse momentum distri- 
vantages: it is a smooth procedure bution of W + W~ pairs in pp collisions 
without cutoffs; the predictions are calculated at NLO, with Herwig (mul- 
guaranteed to be correct at NLO and tiplied a if -factor) and with MC@NLO 
it does not break the resummation of [54]. 

large logarithms. From a practical point of view it has the (minor) draw- 
back of some events with negative weights, however the fraction of negative 
weight events is low (about 10% in the example shown above) and they are 
uniform negative weights, so they should have little effect on the conver- 
gence of the results. Another limitation is that to implement this method 
it is necessary that one understand the Monte Carlo event generator suffi- 
ciently well as to be able to derive the function M(x), i.e. the effective NLO 
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correction already embodied in the event generator. This however is almost 
certainly inevitable: there is no way of ensuring a truly NLO result without 
taking into account what is already included in the event generator. 

6. Conclusions: testing QCD? 

An apology is perhaps due at this stage to those readers who would have 
preferred a detailed discussion of the evidence from final-state measurements 
in favour of (or against) QCD as the theory of hadronic physics. I rather 
took the liberty of reinterpreting the title as 'Tests and perspectives of our 
understanding of QCD through final-state measurements.' Such tests are 
vital if we are to extend the domain of confidence of our predictions, as has 
been discussed in the cases of diffraction and power corrections. 

The tests of course should be well thought through: some considerations 
that come out of the still to be fully understood 6-excess story are (a) 
the importance (as ever) of quoting results at hadron level, not some ill- 
defined parton level; and (b) that if carrying out a test at a given level 
of precision {e.g. NLO), it is necessary that all stages of the theoretical 
calculation (including for example the determination of the fragmentation 
function), be carried out at that same level of precision. 

Another, general, consideration is the need for the Monte Carlo models 
to be reliable and accurate, whether they be used to reconstruct data or to 
estimate backgrounds. This is especially relevant in cases where the actual 
measurements are limited to corners of phase space or where large extrap- 
olations are needed. In this context the recent advances in the extension 
of Monte Carlo models to NLO accuracy is a significant development, and 
in the medium term we should expect progress from the current 'proof-of- 
concept' implementations to a widespread availability of NLO-merged event 
generators. 

To conclude, it could well be that a few years from now, many of the 
measurements and theoretical approaches discussed here will have made it 
to textbooks as 'standard' QCD. We look forward to future speakers on 
this topic have an equally varied (but different) range of 'until recently 
controversial' tests of QCD to discuss! 
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